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ABSTRACT 

Using observed GALEX far-ultraviolet (FUV) fluxes and VLA images of the 21-cm HI column densities, along with estimates of the 
local dust abundances, we measure the volume densities of a sample of actively star-forming giant molecular clouds (GMCs) in the 
nearby spiral galaxy M83 on a typical resolution scale of 170 pc. 

Our approach is based on an equilibrium model for the cycle of molecular hydrogen formation on dust grains and photodissociation 
under the influence of the FUV radiation on the cloud surfaces of GMCs. 

We find a range of total volume densities on the surface of GMCs in M83, namely 0. 1 - 400 cm"^ inside R25 ,0.5-50 cm"^ outside i?25 ■ 
Our data include a number of GMCs in the HI ring surrounding this galaxy. Finally, we discuss the effects of observational selection, 
which may bias our results. 
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• This paper aims to measure the total gas densities in a sample of 

' giant molecular clouds (GMCs) in the nearby spira l M83 (NGC 
5236), using a method initially proposed by Allen e t al.l (Il997h . 
This method is based on the simple (but unavoidable) fact that 
ultraviolet photons from young, newly-formed stars will react 

. back on the surrounding parent GMCs, dissociating the molecu- 
lar gas on their surfaces and turning the (virtually invisible) H2 
into its easily-detected atomic form. The motivation for this ap- 
proach was provided by the discovery of lAllen et alj (1198 6) that 
the HI delineating the spiral arms in a nearby galaxy showed 
a large-scale morphology that was more consistent with pho- 
todissociation near the HII regions than it was with compression 
of the HI farther upstream in the spiral shock. The presence of 
atomic hydrogen is then indicative of a photodissociation region 
(PDR). The metho d was first applied in some detai l to MlOl by 

' ISmith etall (I2OOOI) and more recendy to M81 bv iHeiner et alj 
(I2008h 

The large, nearby spiral M83 has been a frequent target 
of searches for molecular gas using the CO(l-O) spectral line, 
which is easily detectable in M83. However, such studies have 
only recently begun to be carried out with sufficient linear (spa- 
tial) resolution to discern differences in the location of molec- 
ular gas and the emissions from young, hot stars. For example, 
iLord & KennevI (1199 lb describe the presence of molecular (CO) 
emission 300 pc downstream from M83's eastern spiral arm 
dust lane, detected with a 5" x 12" beam. iRand et alj d 19991) 
find that the CO emission is spatially separated from the dust 
lane as well as the young stars on scales of a few hundred par- 
sec. They suggest UV heating, cosmic-ray heating or a two- 
component molecular phase to explain this morphology. Their 
observations show features that show similarity to the largest 
GMCs in the Milky Way with masses on the order of several mil- 
lions of solar masses. IWall et al] d 19931) explore the CO content 
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in the central kiloparsec of M83 and infer a low-density com- 
ponent (n/Zi ^ 10"^ - 10'*cm"^) and a warm (above 50 K), high- 
density component (hh-. > 10^ - lO^'cm^^'). Crosthwaite et al] 
(|2002|) study CO(l-O), CO(2-l) and neutral gas in M83, observ- 
ing a strong truncation of the molecular disk at 6' accompanied 
by a warped atomic outer disk. Using the conventional assump- 
tions about how to convert CO surface brightness to H2 column 
density, they conclude that roughly 80% of the total gas mass in 
M83 is H2. The HI 21 cm and CO surface brightness are found to 
be correlated, but with a complex pattern of offsets. They spec- 
ulate that the temperature of the CO gas is > 20K in the nucleus 
and < 7K in the outer disk. Lundgr en et al.l ([2004) give a full 
overview of previous observations and present full CO (1-0) and 
CO (2-1) maps based on thousands of telescope pointings, but 
with a modest spatial resolution of ~ 1 kpc. They conclude that 
the molecular gas spiral arms mostly trace the dust lanes. They 
expect the H2 mass to dominate that of HI within 7.3 kpc of 
the center. At their linear resolution, CO a nd HI emissi o ns are 
correlated strongly within the optical disk. IVogler et al.l (I2005h 
look into various tracers of star formation, including Polycyclic 
Aromatic Hydrocarbon (PAH) emission lines, another potential 
indicator of PDRs. Various combinations of these lines prove to 
be good tracers of star forming regions in M83. They are found 
predominantly in the spiral arms. We studied th e occurrence of 
PAHs near PDRs in M81 (jHeiner et all l2008h previously and 
confirmed that they are found near the star forming regions in 
almost all cases. Similar M83 data was not available at this time. 

iThilker et al.l (l2005l) find and discuss sites of recent star for- 
mation in the extreme outer disk of M83, associated with the 
warped HI disk of M83, which raises questions about star for- 
mation efficiency and modes of star formation. The nature o f 
this outer disk is investigated further in lGil de Paz et al. I (l2007bh . 
who find evidence that individual young stars are responsible 
for the UV emission in the outer regions. Zaritskv & Christleii^ 
(I2OO7I) find that these UV sources are quite common out to 
2 X R25- In this paper we include candidate PDRs outside the 
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Fig. 1. The location of the candidate PDRs are shown with the GALEX FUV image in the background. Not all FUV sources are 
visible on this image, but all sources are sources of FUV radiation. The inner circle signifies a 1 kpc galactocentric radius, the outer 
one signifies Ris- 



main optical disk of M83 in an effort to probe the amount of 
available gas in this environment. 

The suspected high densities of molecular gas in M83, the 
morphology of the gas in the spiral arms, and the evidence of 
recent star formation in the outer regions, make this galaxy an 
excellent target for a study of photodissociated atomic hydrogen 
and the volume densities of M83's GMCs using our method. 

In ^ we explain our approach and the data we used. In ^ 
we present our results. We discuss and summarize these results 
ingi 



2. Method 

In this paragraph we briefly descri be the data and how w e used it. 
The same method that we used in iHeiner et al. I (I2008h has been 
applied here. For every candidate PDR, the far-UV flux, local 
HI column densities, dust-to-gas ratio and separation between 
the UV source and the surrounding HI features were determined 
to calculate the total hydrogen volume density. 



The distance to M83 is taken to b e 4.5 Mpc, as measured 
from Cepheids by iThim et all (|2003!). We used GALEX UV 
data of M83, where the FUV and near-UV (NUV) images were 
matched for astrometric accuracy. The images have an estimated 
angular resolution of 8", which is a bit poorer than the generic 
GALEX resolution, since M83 is on the edge of the GALEX 
field of view due to pointing constraints. We did not use the 
NUV image otherwise. We obtained the 21 cm HI data from 
die THINGS team (Walter et al., 2008, AJ, submitted). The ra- 
dio beam size is 10.4 x 5.6". The data were analyzed using the 
Groningen Image Processing System GIPSY ( Ekers et al.. . 19731 : 
Ivan der Hulst et al ] |1992l:IVogelaar & TerlouwIEoOlh . 



Although spatial completeness was not our aim, we at- 
tempted to select our candidate PDRs to be bright and isolated on 
the GALEX image and to have a uniform spread in galactocen- 
tric radius as well as FUV flux, in a similar way as was done in 
Heiner etal. (2008). First, GIPSY's POINTS routine was used 
to automatically detect FUV sources with a signal-to-noise of 
20 or more. The other input parameters were kept at their de- 
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Gil de Paz et al. 2007 (model fit) 



Gil de Paz et al. 2007 (high metallicily) 
N ' ■ ^' * • - i - . . Brcsolin & Kcnnicutt 2002 



Glide Paz el al, 2007 (dual) 



Fig. 2. Different possible dust-to-gas ratio models are plotted . 
There are three models suggested by iGil de Paz et al. I (l2007bh 
and Gil de Paz (private communication, 2008), namely their high 
metallicity fit (our default one), their fit based on PDR modeling 
and a slope based on the high metallicity fit inside R25 combined 
with a constant (low) metal hcity in the outer regions. W e also 
considered the one based on lBresolin & Kennicuttl(l2002l) . 

fault value. This yielded 86 source^ Using a histogram of the 
distribution of the sources with galactocentric radius, we ran- 
domly removed a small number of sources to even out the distri- 
bution. We also removed sources that m atched detections in the 
HST Guide Stai- Catalog v2.2, see e.g. lLasker et al.l(l200 8). We 
manually added some sources to improve the spread in galacto- 
centric radius, selected by eye. These sources stood out visually, 
but we did not enforce a specific signal-to-noise ratio. Finally, 
these sources were automatically fitted with gaussians. In some 
two or three cases this led to visible confusion from nearby FUV 
sources as can be seen in the location plot (Fig. [TJ. The final 

source count was 76. These sources are listed in Table [T] 

We used the tilted ring model from Rogsta d et al.l (1 19741) to 
correct for a different position angle and inclination at larger 
radii. This deprojection correction changes the actual galacto- 
centric radii of the candidate PDRs as well as the assumed de- 
projected separation between the UV source and its associated 
HI patches. This inf ormation was use d in combination with the 
metallicity data fromlGil de Paz et"aLl (,2007b,) to obtain the dust- 
to-gas ratio (see llssa et al ■I ll990h . We used 



log^Mo = -0.051/? + 0.43, 



(1) 



where 6/60 is the local dust-to-gas ratio scaled to the value in 
the solar neighborhood. W e adopt a solar metallicity of 8.69 
jAUende Prieto et alj I200I . The galactocentric radius R is in 
kiloparsec. In subsequent plots we normahze R with R2i, w hich 
is 7.63 kpc (or « 6') from I Vila-Costas & Edmundsl (1 19921) . We 
assume that this relation covers the full range of galactocentric 
radii we investigated. It becomes clear from the literature that the 
metallicity in the outer regions of M83 is still uncertain and re- 
sults have only recently become available. For this reason we ex- 
plore a number of alternatives. We also considered, among oth- 
ers, Bresolin & Kennicutt (2002). The different alternatives are 
plotted in Fig.[2l We have adopted the most recent high metal- 
licity fit in lGil de Paz et al.l (l2007b) as our preferred model. 



The (local background subtracted) FUV fluxes of the can- 
didate PDRs were determined next. For each such region, the 
distances from the central source to the nearest HI patches was 
determined. An HI patch is defined as a local maximum in the HI 
column density. A background level of 2 x 10^" cm"^ was sub- 
tracted within 0.9/?25 and 5x 10'^ cm"^ outside this radius. These 
values are based on the general background level of HI emis- 
sion surrounding M83. We also measured HI background levels 
within the optical disk of M83 as well as considered minimum 
column density values throughout M83. This is an attempt to 
isolate the HI formed by dissociating radiation. Since the back- 
ground level is close to the sensitivity limit, some measurements 
end up close to zero. The resulting values of n are generally not 
very sensitive to varying the background HI column. R e writin g 
Equation |2] to yield n (see Equation 4 in iHeiner et al.l (|2008[) ) 
shows that n does not change much as long as the subtracted 
background column remains substantially lower than the scaling 
factor of 7.8 x 10"*^ cm"^, as is the case here. We measured dis- 
tances out to 22". Those distances were then used together with 
the FUV flux to calculate the incident flux impinging on each HI 
patch. A 1500A foreground e xtin ction correction of 0.52 mag 
based on' Schlegel eTaD (Il998l) and lGil de Paz et al.1 (l2007ah was 
applied to the FUV fluxes. Two example sources are plotted in 
Fig.|3] 

Finally, all these elements were used to calculate the to- 
tal hydrogen volu me density by inverting Equation A2 from 
lAllen et al] (|2004|) and inc luding a variable dust-to-gas ratio (cf. 
Equation 6 in lAllenir2004l) . The final equation to be inverted is: 



7.8 X lO^o 
Nhi In 

6/60 



1 + 



IO6G0 / 6 

n \6q 



-1/2 



cm 



(2) 



This translates into a typical rms noise in the FUV flux of 1%. 



where Nhi is the (background subtracted) atomic hydrogen col- 
umn density (in cm"^), 6 1 60 is the dust-to-gas ratio scaled to the 
solar neighborhood value. Go is the incide nt flux mea s ured a t 
the HI patch (the same Gp as i ised i n e.g. lAUen et al. I (12004 . 
Appendix B - see iHeiner et al.l (l2008h for a more extended de- 
scription of how this is measured) and n is the total hydrogen 
volume density, where n - nHi + 2nH^ (cm"-'). This gas is mostly 
atomic on the surface of the GMC and mostly molecular deep 
inside the cloud. 



3. Results 

The measured FUV fluxes are presented first, followed by the as- 
sociated HI column densities. The derived incident fluxes Go are 
presented next, followed by the resulting total hydrogen volume 
densities. Finally, we deal with a number of selection effects. 

The locations of our candidate PDRs are shown in Fig.[T]and 
tabulated in Table [T] We attempted to get a selection of sources 
representative in both their galactocentric radius as well as their 
UV flux, as was mentioned previously. The table also includes 
the aperture of each FUV source within which its flux was deter- 
mined, the starting grayscale level in the detailed plots and the 
grayscale level increment, where each subsequent level is the 
previous level multiplied by the increment. While the starting 
grayscale level is usually 1 x 10 '^cm"^, it is different for some 
sources to enhance the grayscale contrast. (See Fig. |3]for two 
sample plots; the full set is available in the electronic edition of 
this Journal). 

The disk of M83 is warped at larger galactocentric radii. This 
influences the deprojected galactocentric radius that we used 
to determine the dust-to-gas ratio. The distance from the FUV 
sources to the associated HI patches is also affected, since we 
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57= 56= 13^36"'55= 37^^ 13^36"'59= 36"'58' 

R.A. (2000.0) R.A. (2000.0) 



Fig. 3. Sources 4 (left panel) and 24 (right panel) are plotted as examples of candidate PDRs, showing HI grayscales and FUV 
contours. The location of the central source is marked with a black star, the locations of the HI patches are marked with white 
crosses. The FUV contours start at 7 x 10^'^ergs cm"^ s"' A"' and each consecutive level is 1.34 times the previous one. The HI 
grayscale levels are constructed similarly, see Table [T] Both sources show a central source surrounded by HI complexes, where the 
distribution is smoother in Source 24. The full set of sources is available in the electronic version of this paper (without the white 
crosses). 

Table 1. Locations and FUV fluxes of candidate PDRs. The full table is available in the electronic edition of this Journal. 



Source No. 


R.A. (2000) 


DEC (2000) 


Radius 


Fpuv 


Aperture 


Grayscale start 


Increment 








(kpc) 


10"'^ ergs cm"^ s"' A"' 




(10" cm-2) 




1 


13 37 01.346 


-29 52 44.72 


1.1 


5.16 


19 


1 


1.35 


2 


13 36 56.700 


-29 51 54.00 


1.3 


2.08 


16 


1 


1.33 


3 


13 37 03.402 


-29 52 51.43 


1.5 


2.20 


14 


1 


1.37 


4 


13 36 55.827 


-29 52 13.58 


1.5 


2.56 


18 


1 


1.30 


5 


13 36 59.158 


-29 53 02.08 


1.6 


6.83 


21 


1 


1.36 


6 


13 37 02.955 


-29 50 48.84 


1.6 


12.24 


14 


1 


1.37 


7 


13 37 04.572 


-29 50 58.09 


1.6 


6.89 


12 


1 


1.33 


8 


13 37 02.445 


-29 53 05.43 


1.7 


1.18 


10 


1 


1.35 


9 


13 37 08.120 


-29 51 02.97 


2.4 


20.81 


20 


1 


1.37 


10 


13 36 54.317 


-29 53 10.26 


2.5 


21.61 


17 


1 


1.35 



use the position angle and declination to deproject the separa- 
tion Phi (see Heiner et al. (2008) for a more detailed description 
of the procedure and projection issues). TablelHshows the valu es 
we adopted from the tilted ring model o fi Rogstadet all (119741) . 

The background subtracted FUV flux of each source is plot- 
ted in Fig. |4] and listed in Table [T] The background UV level 
was determined at the radius listed in the aperture column. A 
range of fluxes can be seen, which seems to be decreasing to- 
wards the outer parts of M83. However, crowding effects are a 
possible cause of this trend. This is indicated in Fig. |4] In ad- 
dition we have made no attempt to obtain a statistically com- 
plete sample of sources. A detection limit is plotted based on the 
noise level of the FUV data and the method we used to select 
our sources. These fluxes have not been corrected for extinction. 
We did not se l ect sou rces close to the central starburst (See e.g. 
iBoissier et al.l (l2005h for a description of the UV radial profiles 
found by GALEX). 

In an attempt to capture as much information as possible 
per candidate FDR, we determined the distance pHi of every 
HI patch to its central UV source (Fig. [Sj. This means that ev- 
ery FUV source has a number of measurements associated with 
it, which appear in the plot as a number of points at the same 
galactocentric radius. At an arbitrary distance of about 480 pc 



from the FUV source (roughly 22") we stopped recording HI 
patches, since at that distance the association of these patches to 
that particular UV source becomes highly unlikely. This cut-off" 
is arbitrary, since it becomes progressively less likely for an HI 
patch to be associated with the central UV source. The cut-off" 
was chosen with the typical size of large scale PDRs in mind (a 
few hundred parsec) and is labeled 'confusion limit' in the plot. 
Large scale PDRs are expected to span a few hundred parsec 
at most and other UV associations will also be encountered. The 
separations are quantized since we used 1" rings. In certain cases 
the candidate PDRs themselves are fairly close to each other. In 
such a case an HI patch can be attributed to two UV sources, but 
without further morphological information no distinction can be 
made. The minimum observable distance is determined by the 
resolution of our data and is about 30 pc. This plot also shows 
the coverage in galactocentric radius, which shows a few gaps 
where no FUV source was found (or selected) using the selection 
method described earlier. The full results are tabulated in Table 
|3] where the different HI patches per FUV source are labeled a, 
b, c, and so on. pni, Nhi, the incident flux Go, the source contrast 
G/Gbg, the total hydrogen volume density «, and the fractional 
errors are listed in this table. 
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0,0 0-2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

R/R25 

Fig. 4. The FUV fluxes of the candidate PDRs are plotted here 
on a logarithmic scale. Outside of R25 the source fluxes seem 
to flatten out. A detection limit has been included, as well as 
a possible crowding limit. The angular resolution of the data is 
approximately 8". 



In rare cases the HI column density is very low, where the 
measured column is close to the value of the subtracted back- 
ground. No columns higher than * 2.5 x 10^' cm"- were found. 
This is potentially a beam smoothing effect, where we only re- 
solve HI patches partially. We conclude this from our work on 
M33 (Heiner et al., in prep.). At higher linear resolution we ob- 
serve higher HI columns, which points to better resolution of the 
HI clouds in M33. 



Table 2. Adopted tilted ring values from lRogstad et alj ( |1974|) . 



Radius (') 


p.a. (degrees) i (degrees) 




0.5 


45 


24 


1.0 


45 


24 


1.5 


45 


24 


2.0 


45 


24 


2.5 


46 


24 


3.0 


46 


24 


3.5 


46 


24 


4.0 


46 


25 


4.5 


46 


26 


5.0 


49 


27 


5.5 


52 


28 


6.0 


59 


30 


6.5 


67 


32 


7.0 


75 


34 


7.5 


83 


36 


8.0 


91 


37 


8.5 


99 


39 


9.0 


106 


40 


9.5 


113 


42 


10.0 


120 


43 


10.5 


127 


45 


11.0 


133 


46 


11.5 


139 


47 


12 


145 


49 


12.5 


150 


50 



• • • 
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o* • • 

•• • • • 

• • 

o • • • • 

• • • • • 

resolution limit 



The incident flux Go (corrected for foreground extinction) 
on each individual HI patch is plotted in Fig. |6l A set of points 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
R/R25 

Fig. 5. The distances pni from the central UV source to individ- 
ual HI patches are plotted here. Each source can have multiple 
HI patches associated with it. As an example, 10 HI patches as- 
sociated with one particular UV source are plotted as open cir- 
cles here, at R/R25 ~ 1 (Sources 61 a-j in Table[3l three sources 
have an identical separation). At a separation larger than 22" 
no HI patches are recorded anymore. In the outer parts fewer 
HI patches are detected, possibly because of sensitivity issues. 
Some gaps in the radial coverage can be seen, where no sources 
were found with our selection method. 



in the vertical direction shows all resulting Go for an individual 
UV source with multiple HI patches. The average Go declines 
going outward, with no distinct break at the edge of M83's opti- 
cal disk. A hint of an upturn at the location of M83's HI ring is 
visible at the extreme right end of the plot. As Go is a combina- 
tion of the FUV flux and the FUV source / HI patch separation, 
it is also influenced by the possible crowding effects that are in- 
dicated in Fig. m The minimum Go (0.005) is due to the maxi- 
mum Phi that we measure, combined with the lowest measured 
flux, and is therefore a sensitivity limit. The highest Go (565) 
that we could obtain is determined by the lowest measurable pHi 
and the highest FUV flux. Both limits are plotted. The open cir- 
cles are incident fluxes with a source contrast of less than 0.5. 
This means that the incident source flux has 50% of the strength 
of the ambient FUV field incident on the HI patch or less. The 
HI corresponding to the closed circles can be assumed to have a 
higher chance of having been produced by photodissociating UV 
radiation, since the central UV source is dominant with respect 
to the background radiation field (the source contrast is higher). 
However, these highest source contrasts occur when the separa- 
tion between the central UV source and the associated HI patch 
is the lowest - at close separation projection effects are more of 
a problem, which means a higher chance of underestimating the 
separation and overestimating the incident flux. A high source 
contrast at an HI patch therefore needs to be viewed with cau- 
tion. 

The (beam averaged) atomic hydrogen column density of 
each HI patch in our candidate PDRs is plotted in Fig. |2l A 
selection of local maxima are counted as HI patches. Near the 
center of M83 there is little HI emission visible against the con- 
tinuum source at the nucleus (e.g. Tilanus & Allen 1993). The 
maximum observed column density rises to a broad peak around 
0.5 R25- The end of M83's disk is clearly visible where the max- 
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maximum of the method 




Ti h« 



R/R25 

Fig. 6. The incident flux Go is plotted here. Some holes in the 
radial coverage are visible. The fluxes within 0.27?25 seem to 
be lower Generally, fluxes are dropping going outward. The 
maximum of the method consists of the Go that would have re- 
sulted from our maximum measured FUV flux, combined with 
the lowest separation that can be resolved at the given resolu- 
tion. Sources with a source contrast of less than 0.5 are plotted 
as open circles. 



O 15 

A 



V 



optically tliick column 



m Km 



sensitivity cul-off - 



0,0 02 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 

R/R25 

Fig. 7. The (beam averaged) HI column densities are plotted 
here. Values range from around the sensitivity limit at 1 x 
10^** cm"^ to almost 2.5 x 10^' cm"^. Columns are highest within 
R25 and a small upturn can be seen signifying M83's HI ring at 
about 1.7 /?25- 



imum HI column densities drop, but HI is still detected. Finally, 
there is a slight increase where M83's HI ring is located. 

An HI background value was subtracted as described previ- 
ously, in an attempt to capture only the HI produced by the inci- 
dent FUV radiation from the FUV source nearby. The resulting 
values of n are relatively insensitive to varying the background 
HI column. Background-subtracted values below the sensitivity 
limit were discarded. 

The calculated total hydrogen volume density is shown in 
Fig. [8] based on the Gil de Paz et al. (2007b) high metallicity as- 
sumption, which is our preferred dust model. Values span a range 
of two orders of magnitude mostly, with values of a few thousand 




0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 
R/R25 

Fig. 8. The total hydrogen volume density is plotted, assuming 
relatively high metallicity (our preferred dust model, Equation[T]) 
outside the optical disk. Closed circles indicate a source contrast 
higher than 0.5, open circles indicate a lower source contrast. An 
approximate range of values of n is indicated (dashed lines): 0.1 
- 400 cm"^ inside R25, 0.5 - 50 cm"^ outside R25. 



cnT^ in the inner parts of M83. Outside M83's optical disk the 
values of n are generally lower In this plot, again closed circles 
are plotted for HI patches that show a source contrast higher than 
0.5 (meaning the source is half as bright as the UV background 
or stronger), and open circles for those with lower source con- 
trast. The former may be considered more reliable than the latter, 
since the UV source is considered to be more clearly responsi- 
ble for the dissociated HI. The FUV source contrasts G/Gtg are 
listed in Tabled The frac tional errors in Table[3]are as calculated 
in the lHeiner et al. I (l2008l) M81 results, except that we used a rel- 
ative error for the dust-to-gas ratio of 20% instead of a constant 
error The average fractional error is 0.7. Since photodissocia- 
tion rates are highly sensitive to the local dust-to-gas ratios, we 
looked at a number of recent results for M83. 

Different alternative total hydrogen v olume densities based 
on va riou s dust models s uggested by Bresolin & Ke nnicuttI 
(I2OO2I) and lGil de Paz et all (|2007b) are shown in Fig.|9] The dif- 
ferent du st models can be found in Fig .[2l The metallicity model 
based on lBresolin & KennicuttI (|2002|) uses 



log^Mo = -0.0217? -H 0.17. 



(3) 



The slope of this model is relatively shallow. The dual metal- 
licity model in the middle panel uses Equation [T] except that 
6/60 = 0.15 for all PDRs outside of R25, assuming low and con- 
stant metallicities in the outer regions of M83 (Gil de Paz, private 
communication, 2008). The bottom panel uses 



log6/6o = -0.1007? -t- 0.85. 



(4) 



The slope of the metallicity is relatively steep in this case since a 
single slope is assumed, with high metallicities inside the optical 
disk of M83 and low metallicities outside. 

The top panel densities are similar to the ones in Fig. [8] since 
the underlying dust models are also very similar. The dual metal- 
licity model in the middle panel yields densities in the outer re- 
gions that appear to be unusually high. Finally, the densities in 
the bottom panel are extremely depressed at small galactocen- 
tric radii, a result of the high dust content. Insufficient data are 
available to pick one dust model over the other We continue 
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Fig. 9. The total hydrogen vol ume density is plot t ed for 
three dust models de ri ved fro m iBresolin & Kennicutt! (l2002h 
and 'Gil de Paz et al.' ('2007b') . The results based on the 
[Bresolin & Kennicutt (2002) dust model are similar to the one 
from lGil de Paz et al.l (|2007b) assuming a high metallicity out- 
side the optical disk (Fig.lS). In the dual metallicity case, the hy- 
drogen densities are relatively high in t he outer regions. The sin- 
gle slope fit based on PDR modeling bv lGil de Paz et al.1 (l2007bh 
is steep because of relatively high metalhcities near the center 
of M83 and low metalhcities in the outer regions. Densities are 
markedly lower in the center in this case. 



4.0 
3.0 

2.0 



0.0 
-1.0 



0-0-5 0.5-1 1-1.4 1.4-1.8 

R/R25 bin 

Fig. 10. Box-and-whisker plots for the total hydrogen volume 
densities in Fig. [8] are displayed, where the data has been di- 
vided into four galactocentric radius bins. The boxes span the 
lower to upper quartiles of the data (50% of the datapoints in 
the bin combined), with a line indicating the median. The so- 
called 'whiskers' show the range of the data (1.5 times the inner 
quartile range of points), with outliers plotted individually {+). 
The approximate ranges in Fig. [8] are in good agreement with 
this. The maximum densities in the outer regions are lower than 
the ones in the inner regions of M83 and the range of values is 
narrower 



our a nalysis using the high metallicity model of iGil de Paz et al.l 
(l2007 b). as shown in Fig.|8] 

In order to further quantify the volume densities in Fig. [8] so- 
called box-and-whisker plots are shown in Fig. [10] Since we did 
not aim to select a complete sample of candidate PDRs, we are 
mostly interested in the range of volume densities. The galacto- 
centric radius bins in this plot give a good impression of the dis- 
tribution of the volume densities. There is a wider spread of val- 
ues in the inner regions of M83, accompanied by higher density 
values. The values in the outer regions are more concentrated 
and the maximum densities are lower. This is mostly due to the 
narrower spread of HI column densities outside /?25. Selection 
effects also limit the range of observed densities. 

Our results are summarized in Fig. [TT] which is suitable to 
show the various limiting selection effects. The candidate PDRs 
at a galactocentric radius larger than R25 are plotted as open 
squares. The Roman numerals I through V indicate the differ- 
ent limits of our observations: 

I. The HI column density upper limit of 5 x 10^' cm"^ at a 
typical S/S() of 1.8, due to the HI column becoming optically 
thick. However, the actual upper limit of the HI is closer to 
2.5 X 10^' cm"^ because of beam smoothing effects. These 
values occur around a typical 6/60 of 2.4 and are delineated 
by an extra dashed-dotted line. 

II. The HI lower hmit of 1 x 10^" cm"^ at a typical 6/6q of 1. 
This is a sensitivity limit. 

III. The lowest volume density that we can observe as a con- 
sequence of the radio beam diameter and the HI sensitiv- 
ity limit. For M83, it is approximately 0.2 cm"-'. Points may 
end up lower because of high metalhcities. Ideally, this limit 
would contain the GMC cloud size. Since this is unknown, 
the beam diameter is a reasonable substitute. 




-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 



log Go 

Fig. 11. Our results can be summarized in a plot of Go against n. 
The numerals indicate observational selection effects discussed 
in the text. In this plot, the open squares indicate regions outside 
ofRis- 

IV. The minimum Go (0.005) that we can use (at a source con- 
trast as low as 1%). 
V. The maximum Gq (565) that we can obtain from the data. 



4. Discussion and Conclusions 

In this paragraph we compare our results to H2 densities based 
on CO measurements. Then we discuss background issues, ex- 
tinction and dust issues. We end with a summary of our findings. 
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Table 3. HI Measurements, Incident Fluxes and Volume Densities 



Source No. 


Phi (pc) 


Nhi (10"' cm--) 


Go 


G/Gbg 


n (cm ^) 


Fractional error 


la 


196 


1.26 


1.41 


2.05 


120 


0.8 


lb 


349 


2.77 


0.45 


0.65 


12 


0.7 


Ic 


371 


4.28 


0.40 


0.57 


5 


0.7 


Id 


436 


6.99 


0.29 


0.41 


1 


0.9 


le 


458 


6.69 


0.26 


0.38 


1 


0.9 


If 


458 


9.40 


0.26 


0.38 


0.3 


1.2 


2a 


262 


1.41 


0.32 


0.57 


24 


0.7 


2b 


262 


2.62 


0.32 


0.57 


10 


0.7 


2c 


284 


2.77 


0.27 


0.49 


8 


0.7 


2d 


371 


2.47 


0.16 


0.28 


5 


0.7 


2e 


393 


2.62 


0.14 


0.25 


4 


0.7 


2f 


393 


2.92 


0.14 


0.25 


4 


0.7 


2g 


415 


3.67 


0.13 


0.23 


2 


0.7 



4. 1 . CO results comparison 

The study of iRand et al.l (Il999h is detailed enough to find typi- 
cal sizes and masses of GMCs in the eastern arm of M83. These 
results can be compared to ours with some additional assump- 
tions. Using a CO(l-O) map with kinematic information, they 
find CO masses and virial masses onindicatingindicating the or- 
der of 1 X 10^ Mq. For example, their source 9 (their Table 4) 
yields a mass derived from CO emission of 7.3 xlO''Mo and 
3.6 xIQ^Mq for the virial masses, at a 50-80% uncertainty. This 
source corresponds to our FUV source no. 13. Since our method 
does not yield any information on cloud sizes (we only observe 
HI on the surface of GMCs), nor on kinetic temperatures in these 
clouds, we will take typical values of our measured volume den- 
sities and compare densities using a typical GMC radius of 75 
pc for a spherical cloud, n = 2n/fj inside the GMC and a typical 
density near our source no. 13 is 20 cm^^, or = 10 cm"^. 7.3 
xlO^Mo is equivalent to about = 85cm-"', which is in rea- 
sonable agreement with our results considering the uncertainties 
in both results (70-80%). More recent unpublished data (Lord, 
private communication, 2008) may be quantitatively different 
and more detailed, but does not yield substantially different re- 
sults at the resolution of the data that were used in this paper 

One of the key differences between M81 and M83 is the 
abundance of CO emissio n. While CO is extr emely faint and 
hard to detect in M81 (e.g. Casasola et"aLll2007 ), M83 displays 
bright CO features. The method we used here to find hydrogen 
densities (and molecular hydrogen densities in GMCs) does not 
yield any morphological information, since we are observing HI 
on the surface of PDRs. We can note, however, that the size and 
scale we assume for our candidate PDRs (the range of values 
of Phi) is consistent with findings using CO. Since we did not 
find suitable candidate PDRs in the inner kiloparsec of M83, 
we cannot compare our results with for example the results of 
IWall et al] ([T993), although the densities we find match their 
low-density component. A s to the truncation of the molecular 
disk that was mentioned by Crosthwaite et al.l (l2002i) . we do see 
a difference between the inner and outer disk of M83, but it does 
not seem to be that significant in our results. No differences are 
seen in candidate PDRs in arm or inter-arm regions, insofar as 
the large scatter in our results allows us to draw this conclusion. 
We do not find m ore than a handful of regions with the den- 
sities indicated by iHosokawa & Inutsukal (l2007l) . although the 
expected sizes of our candidate PDRs are the same. Finally, the 
hydrogen densities presented here are similar to the M8 1 results 
but extend to higher maximum values. This is consistent with 
the brighter CO emission in M83, since the higher H2 densities 



will lead to a greater degree of excitation of the CO molecules. 
Alternatively, this could also i ndicate a higher fract ion of CO 
molecules in M83's clouds (e.g. lCasasola et af]|2007h . 

4.2. Background levels, extinction and dust 

The FUV fluxes seem to be higher in the inner parts of M83 and 
decreasing (on average) going outward, as is shown in Fig. |4] 
These higher fluxes in the inner parts could be caused by brighter 
sources, and/or a larger number of sources. Since the individ- 
ual sources are unresolved, no firm conclusions can be drawn 
here. Outside R25 the measured fluxes are roughly constant. At 
the same time, the abundances of PDR-produced HI surrounding 
these FUV sources (Fig.|7]- note that these are individual mea- 
surements, not annular averages) seem to follow a similar trend, 
consistent with its connection to the photodissociating UV radi- 
ation. The FUV source contrast does not vary with galactocen- 
tric radius, so while both source and background radiation field 
decrease in intensity towards the outer regions of M83, their rel- 
ative strength does not change. 

Another factor that could influence the FUV radiation is 
galactic foreground extinction. The higher the applied fore- 
ground extinction, the higher the resul ting total hydrogen vol- 
ume density. The 0.52 mag extinction jSchlegel et aljfl998h to- 
wards M83 is significantly less tha n the 1.37 mag we used to- 
wards M81 in iHeiner et al.l (l2008h . The ISchlegel et al] (Il998l) 
extinction correction for M81 would be 0.58, significantly low- 
ering the total hydrogen volume densities. We find a wider range 
of values in M83, including higher gas densities. We therefore 
expect M83 to harbor more gas than M8 1 . The M8 1 and M83 re- 
sults are hard to compare directly because the extinction correc- 
tions were based on different sources in the literature. We intend 
to compare the two galaxies, together with similar M33 results, 
more quantitatively in a future paper, with consistent extinction 
corrections. 

Our results are also sensitive to the dust-to-gas ratio. The 
slope of M83's metallicity is relatively shallow. Our preferred 
dust model assumes that metallicities in the outer parts of M83 
remain high. Lower metallicities would result in higher gas den- 
sities. Our results in these regions are therefore most likely lower 
limi ts and much more g as could be present. The recent results 
by Gil de Paz et al.l (l2007bi) are ambiguous in the sense that the 
authors provide high metallicity and low metallicity results, de- 
pending on the adopted model. In the high metallicity scenario, 
the slope of M83's metallicity remains shallow out to large 
galactocentric radius. The low metallicity scenario is accompa- 
nied by a sharp drop in metallicity starting at approximately R25 . 
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ISmith et al] (l2000h found a similar range of gas densities in 
MlOl using basically the same method with additional extinc- 
tion corrections. Their candidate PDRs were all within MlOl's 
R25, but the FUV luminosities within that range are similar As in 
our previous M8 1 results, no internal extinction correction was 
applied. If any such correction were applied, it should scale with 
the dust-to-gas ratio. This means that the extinction would de- 
crease going outward. The FUV fluxes would be even higher in 
the center of M83 (leading to higher total hydrogen volume den- 
sities), while the fluxes in the outer regions of M83 would not be 
affected much. 

4.3. Summary 

In summary, we have investigated atomic hydrogen found in can- 
didate PDRs in M83 and used the physics of these PDRs to de- 
rive total hydrogen volume densities. We carefully considered 
the contribution of observational selection eff'ects to our results. 
We find a range of densities: 0.1 - 400 cm inside R25, 0.5 - 50 
cm""* outside R25, based on measurements of HI believed to be 
produced in large scale PDRs in M83. The higher GMC volume 
densities which we find within R25 correlate well with the pres- 
ence of bright CO(l-O) emission there. This points to enhanced 
collisional excitation as one reason for the CO emission in these 
GMCs. We note that this is also consistent with our results in 
M8 1 , where we find little evidence for high GMC volume den- 
sities, and for which the CO(l-O) emission is faint. 

Our measurements go out to a galactocentric radius of 13 
kpc (deprojec t ed). O ur study used the tilted ring model from 
iRogstad et al.l (1 19741) to get proper galactocentric radii outside 
the M83 optical disk. Our results are notably sensitive to the lo- 
cal dust-to-gas ratio, especially since the metallicity in the outer 
regions remains uncertain. 
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